Bosonization of Vector and Axial- Vector 
Currents in 3 + 1 Dimensions. 

CD. Fosco* 
Centro Atomico Bariloche 
84 00 Bariloche, Argentina 

F.A. Schaposnik^ 

Departamento de Fisica, Universidad Nacional de La Plata 
C.C. 67, (1900) La Plata, Argentina 

February 7, 2008 



Abstract 

The bosonization of a massless fermionic field coupled to both vec- 
tor and axial-vector external sources is developed, following a path- 
integral approach. The resulting bosonized theory contains two an- 
tisymmetric tensor fields whose actions consist of non-local Kalb- 
Ramond-like terms plus interactions. Exact bosonization rules that 
take the axial anomaly for the axial current into account are derived, 
and an approximated bosonized action is constructed. 
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Any physical system is suitable of many different mathematical descrip- 
tions, i.e., there is a lot of freedom in the choice of variables used to define its 
configurations. Different choices of variables are equivalent in the sense that 
they describe the same system. An extreme manifestation of this appears 
in some two-dimensional models, which can be described in terms of either 
fermionic or bosonic variables. The equivalence between these two formu- 
lations is made explicit by the so called bosonization rules, that map the 
fermionic to the bosonic variables. 

There has been some progress in the program of extending, at least in an 
approximated way, the bosonization procedure to theories in more than two 
dimensions [0-0. In this letter we are concerned with the problem of finding 
the bosonization rules for a simple system in four space-time dimensions, that 
of a massless fermionic field in the presence of vector and axial sources. This 
will allow us to find the bosonization rules for the vector and axial fermionic 
currents, this last one being consistent with the axial anomaly [ID in the 
presence of an external vector gauge field. 

Our starting point is the generating functional for a massless Dirac field in 
3 + 1 (Euclidean) dimensions, coupled to Abelian vector (s^) and axial- vector 
(tfj,) external sources 

Z{s^t^) = J VipVip exp -5($,V>;s M ,^ 

5(^,^5 = -i J ' d 4 x - f -75$^ (!) 

where we have adopted for the 7-matrices the following conventions: 

7j = 7 M , 7s = 7s , {7 M , 7i/ } = 2 V • ( 2 ) 

The addition of the source is motivated by the reason that, in four 
dimensions, the vector and axial currents are independent fermionic bilinears. 
Thus not all the information provided by Z{s^ t^) can be obtained from, say, 
Z(Sfj_, 0). In two dimensions, because of the smaller number of generators for 
the Dirac algebra, these two current are related, and the bosonization rule 
for one of the currents also implies the proper rule for the other. 

By a natural extension of the procedure followed to obtain the bosonized 
version of the generating functional with just a vector source (see for example 
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[0]"@ f° r a detailed description of the approach) , we perform in ([I]) the 
change of variables 

tjj(x) = e^^^VO) , tp(x) = ip'(x)e- ieix) - i ^ aix) . (3) 

In terms of the new variables, the generating functional (0) reads 

Z{s tl ,t ll ) = f V^Vip J(a; s^tp) exp -5(^,^5 <V + ^0, ^ + 9 M a) (4) 

where the primed fermionic fields have been renamed as unprimed for the sake 
of simplicity, and J is the anomalous Jacobian corresponding to this fermionic 
change of variables, a well-known consequence of the chiral anomaly [[11 



This Jacobian is evaluated by using the standard Fujikawa's recipe. To de- 
cide whether the consistent or covariant regularization is to be used in the 
calculation of this determinant, it is important to recall (see for example 
ref.|12[) that the consistent regularization, besides assuring the validity of 
Wess-Zumino consistency conditions, automatically enforces the conserva- 
tion of the vector current. The covariant regularization, in turn, guarantees 
a gauge- covariant form for the anomaly, at the expense of generating an 
anomaly for the vector current as well. Regarding the case at hand, it is 
desirable, although not mandatory, to have a regularization that assures the 
conservation of the vector current, since one has in mind applications of the 
bosonization recipe to situations where interactions involve the vector cur- 
rent, which should then be non- anomalous, while the axial current can have 
an anomalous divergence without spoiling the consistency of the theory. This 
sets the natural choice of regularization in this case to be the consistent one, 
which we adopt. This justifies a posteriori the fact that we have written J in 
(fD as a function of a only, rather than depending also on 6. To obtain the 
Jacobian for this finite transformation, we shall use the techniques described 



in ref. |13| . The consistent regularization of the anomalous Jacobian may be 
based on the hermitean operator 

V = i - ft + 2 75 f (5) 

obtained by performing the analytic continuation —>■ —it^ in the operator 
appearing in the kinetic term of the fermionic action. One then introduces a 
parameter r e [0, 1], to define the r-dependent transformation 

${x) = tfj' r {x)e ird{x) - r ^ a{x) (6) 
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which, of course, induces a change in the operator T> 



-p _ e ir9(x)-r"/ 5 a(x)<j~) e -ir6(x)-ra(x)"/5 

= i - (/-r 09) + Hs(1/-r 0a) . (7) 

The Jacobian J for the finite transformation @ corresponds to r = 1, and 
may be obtained by multiplying the Jacobians corresponding to infinitesimal 
steps dr, each one calculated with the appropriated regularization operator 
V r 

J = exp(2 / dr lim Tr[aj 5 e~ w ?]\ . (8) 

The factor 2 in the exponent comes from multiplying the two Jacobians, 
associated to the changes in the measures T>ip and T>tjj. Although each one 
of these Jacobians does depend on 9, like terms do cancel in the product. 
This is a consequence of the consistent regularization, which uses the same 
regularization operator for both Dip and Dip. 

There just remains to insert into (|S[) the result for the functional trace, 
which can be read from the formula for the anomaly for an infinitesimal 
transformation 



lim Tr 



-tv 2 



1 1 

■ ^A t ^p°"( dnS u 3pS<j -\- — dpt u Opt a ) (9) 



which, in this Abelian case, becomes r- independent, rendering the integra- 
tion over r in (^) trivial. After undoing the analytic continuation on t, the 
anomalous Jacobian for the transformation (H) becomes 



J(a;Sp,tp) = exp 



4vr 2 



J d x cx{x) €pyp a {dpS v dpS a -\- ^dpt v dpt a/ 



(10) 



The next step in the bosonization procedure follows from realizing that, as 
@ is a change of variables, Z cannot depend on either 9 or a. Thus 9 and 
a can be integrated out without other effect than the introduction of an 
irrelevant constant factor in Z, which we ignore. Integration over 9 and a is 
equivalent to integration over two flat Abelian vector fields 9p and a p. 

0^ = dp9 , ap = dpa 
Uv(6) = dp9 v -d v 9p = , fp U (a) = dpa u -d u a>p = 0. (11) 
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The fermionic action obviously depends on a and 9 only through a,p and 9p 
defined in ([0]). By an integration by parts, one sees that so does the Jacobian 
([TP]) • Had the situation been different (i.e., through the dependence of J on 
a and not only on its gradient a M ), it would not be possible to pass from a 
description in terms of a and 9 to one in terms of otp and 9 p. To make this 
explicit, it is convenient to define 
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J d X Oip(x) £ pui pa[s u d pS a + ^t u dpt a j 



12) 



J(ap]s^tp) = exp 
Thus (||) becomes 

Z = J V6rVa»TtyTty 8\J^(6)]8\J^{a)\ 

exp -S{^)^\ s M + 6^,^ + a M ) . (13) 
Formally integrating out the fermionic fields and making the shift of variables 

9^ ^ 9p - Sp , ap — > otp - tp , (14) 

(HD leads to 

Z(sp,tp) = J T>9pVap5{fp U (9 - s)}5{fp V (a - t)] J(ap - tp; Sp,tp) 

x det O + * ^+«7 5 ^) • (15) 

By analogy with the bosonization procedure followed in 0-0, we exponenti- 
ate the two functional delta functions in fll~5|) using two antisymmetric tensor 
fields Ap U and Bp U as Lagrange multipliers 

Z{sp,tp) = J VAp V VBp V V9pVap J(ap - t M ; s M ,t M ) 

exp U y d 4 x[e Mi/pCT A^(9 p 6' CT - 9 p s CT ) + tp Vpa B^ v (d p a (7 - fyO] j 

x det{$ + i $ + . (16) 

The bosonized form of Z can then be obtained by integrating out 9 p and a p 
in (|i~6|) . This produces a generating functional with the tensor fields Ap V and 
Bp V as dynamical variables. This step requires the evaluation of the fermionic 
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determinant, which in four dimensions is necessarily non-exact. Before em- 
barking on such calculation, we derive the rules that map the vector and 
axial- vector currents into functions of the bosonic fields A^ u and B^ u . This 
correspondence requires no approximation and may well be called 'exact'. 
These rules follow from elementary functional differentiation 

3n = $ltfl>) = log^| Sfl =o = -epvpodvAp,, (17) 

5 i 
hti = (V'TsT/^) = -ijjT lo g^l<M=o = -^iivpad v B pa - — e^paSydpSv . (18) 

From the antisymmetry of the tensors A^ v and B^ u , we are entitled to derive 
the equations for the divergencies of the currents: 



dpjp = 



9 A = -^ F M F M- ( 19 ) 



with Fp V = {l/2)e,j, vo ,pF a p. We then see that the bosonization rule (|TJ) 
correctly reproduces the axial anomaly. 

As stated above, although the bosonization recipe dl7|)-(|l8|) for associat- 
ing the fermionic currents with expressions written in terms of bosonic fields 
is exact, the bosonic action governing the boson field dynamics cannot be 
evaluated in an exact form in d > 2 dimensions. Different approximations 
for computing the fermionic determinant would yield alternative effective 
bosonic actions valid in different regimes. We shall describe here the evalua- 
tion of the fermionic determinant in (|16|) to second order in the fields 8^ and 
a p. The use of this quadratic approximation can be motivated by the same 



kind of arguments (see in particular the 'quasi-theorem') used in ref. [14 . 

As usual, it is convenient to work in terms of W, the generating functional 
of connected Green's functions of the fermionic current 

det(0 + i? + ij 5 $ = exp [W{9^ap)} . (20) 

The unregularized form of the quadratic part of W in fl20|) becomes 

Wie^a,) ~ ^Tr [r^ + 75 f) r\9+l* fi] • (21) 



6 



W{6^ atfj) will consist of three parts 



W(9^,a, 



We 



'Hi 



Qp) + WaaioCn, OCfi) + We a (Q^ a t- 



(22) 



corresponding to the three terms in the quadratic part of W. Before ex- 
tracting the three parts of W, a regularization should be introduced, since, 
at it will become clear next, it does have a non-trivial effect. Our choice of 
a regularization is dictated by the requisite of compatibility with the regu- 
larization of the anomalous Jacobian, which we decided to be consistent. A 
form of assuring consistency is to use Pauli-Villars regularization in a way 
that treats the vector and axial-vector vertices symmetrically. This amounts 
to defining the regulated W by 



W reg (9 p ,aJ = -£C s Tr + M s ) -1 (^ + 75 

1 3 = 



+ M S 



where 



M 2 



Cq — 1 G\ — 1 C2 
= /1 2 M\ = 2A 2 - /j, 2 



Ml 



A 



+ 75 ft 
(23) 



(24) 



where A is the cutoff and we have introduced an IR regulator fi that gives 
the gauge field a small mass inside the loop. The regulated version of (f22|) 
becomes 



W reg (e„a^ = W r e °J6„6,) + W r r g (a„a,) + W^J9„a 



(25) 



with 



W aa (a a 



-£C s Tr \(@ + M s )- 1 W + M,)- 1 ? 



s=0 
2 



2 S C ^ 

Z s=0 



+ M S 



M, 



s-l 



(26) 



(27) 



for the first two terms. Knowledge of the mixed term involving 8 P and 
is not necessary, since because of parity violation and Lorentz invariance its 
form is restricted to be 

W r e : g (8„ a,) = jd 4 x J d A y e^ pa d,d u (x)H(x - y)d p a a (y) (28) 
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which is a total derivative, and can thus be safely ignored (at least in this 
Abelian case we are dealing with). We then consider the terms ( p6|) and 
(|2"7D. The first one is evidently equivalent to the quadratic part of the ef- 
fective action for a Pauli-Villars regulated fermionic field in the presence of 
an external vector field 9^. Setting the renormalization conditions at zero 
momentum, we have for this renormalized two-point function 



Wre 9 (^,8») = o / d 4 xd 4 y9,(x)S^F(x-y)9 u (y) 



where 8j^ v is the transverse 5 



J fJ,U 



(29) 



(30) 



and 



F(k) 



F(x - y) 
k 2 



d 4 k 



2tt 2 



(2vr) 4 
dxx{l — x) log 



e iHx ~ y) F{k) 



l + x(l 



XI 



(31) 



As it will become evident at the end, to impose a renormalization condi- 
tion on W ee is equivalent to set the same renormalization condition on the 
(one-particle irreducible) current-current correlation function. Of course, the 
current is the vector one for 9^ and the axial- vector one for a. 



Regarding the term fl27|) , we may rearrange it in such a way that it shows 
a part which is functionally identical to Q2l|), but with as argument, plus 



a remainder AW™ 



ctot 
reg 



where 



s=2 



AW aa 

" reg 



E C * M s Tr 



ft" 1 



(32) 



(33) 



s=0 



-d 2 + M 2 )- 1 fi{-8 2 + M 2 

Of course, the renormalized term W ed (a fJi ,a fJi ) will be given by 

W B \a^ O = - J d 4 xd 4 y a^(x)5^F(x - y)a v (y) (34) 
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with F as in (0). 

A straightforward evaluation shows that AW™ is given by 

= ~7^~v? I d 4 xa^x)a^x) . (35) 



The meaning of this quadratic divergence is that the consistent regularization 
violates axial gauge invariance, as expected. The renormalization of this 
divergence requires the introduction of a mass counterterm, and also the 
choice of renormalization conditions for W aa which should say which is the 
value of the renormalized mass for a. Thus the renormalized AW aa will 
correspond to a finite mass term: 

,2 



in f 

AW M = - — J d 4 xa,(x)a,(x) 



4 ; d * x 



a:, 



x)8L a n( X ) + a ni X ) S L a fi( X ) ( 36 ) 



where 5^ = d ll d~~ 2 d u , and m is the renormalized mass. Of course we can set 
the value of the renormalized mass to zero to this order. If a M is absolutely 
non-dynamical, it is possible to set the renormalized mass equal to zero 
to all orders, since the only correction to the mass term is produced by 
the term we are considering. If the field were dynamical, the anomalous 
behaviour of the axial symmetry would have spoiled the masslessness of a 
to higher orders. From the previous results it follows that we can write the 
renormalized functional W as follows: 

W(9„a,) = ~ J d A xd A y [9,{x)6^F{x - y)6 v (y) 
+ ^(x)5^G(x - y)a v (y) + m 2 a^(x)5l u 5(x - y)a v (y)] (37) 

where 

G(x - y) = F(x - y) + m 2 5(x - y) . (38) 

Inserting Q57p into flTE|), we see that the functional integral is Gaussian 
with respect to 9^ and ot^\ 



Z(s^t^) = JvA tlu VB tliy V9 fl Va tl 
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exp 
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exp J d 4 x[e p>vpa A IJiV (d p 9 CT - d p s a ) + e pupa B pu (d p a a - d p t a )]\ 



exp 



d*xd*y [e^x)6±,F(x-y)e v (y) 



+ a fl {x)5^G{x -y)a u {y) + m 2 a p {x)5lj{x - y)a h 
Performing the Gaussian integration over 9 p and a p in ( |39| ) we get 

Z(sp,t p ) = exp[C(s M ,£ M )] J VA pu VB pu x 

exp ^ i J d x\s p t pv piJ d v A pa -\- t p (a pvpcT 3 u B pcT -\- ~^^a pup(T s u dpS < j^~^ 
d 4 xd i y[A pvp (x)F' l {x - y)A pup (y) + 



exp 



B pup (x)G 1 (x - y)B lw/> (y)]} x exp / d 4 xd 4 yd p B up (x)x 

G" 1 (x - y)S IM/ p ta j3 1 (s a dps^ + ^tadptj) 



where 



.4 
5 



Jpup,a/3f 



d p A up + d v A pp + SpA^ 

pi' 



det 



^ <Vy 
^i/a ^i/jS ^1/7 

\ O pQ Op/3 5 P1 



and 



2(2tt) 



d xd y \ [s p (x)d l/ s x (x) + -t p (x)d u tx(x)] 



Sp,up,aPjG 1 (x - y)[s Q (y)^s 7 (y) + -t a (y)d/3t 1 (y)} 



2(2tt) 
1 



^ | d*xd i yg(x)d- 2 G- 1 (x-y)g{y) 



+ 



2m 2 (27r) 4 

where £ = e pup x(d p s u dpSx + \d p t u dpt\). 



d xd yG(x)d (x - y)Q(y) 
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Summary and Conclusions 



We have applied the bosonization technique developed in ref's [5-9] to the 
case of massless Dirac fermions in four dimensions in the presence of both 
vector and axial-vector sources. This has allowed us to find the bosonization 
rules for both fermionic currents, eqs. (|T7|) -(|T§|), in terms of Kalb-Ramond 
bosonic fields. While the bosonization rule for the vector current can be 
written in a natural and compact form, reminiscent of the well-known two- 
dimensional bosonization rule, 

i>l^ip -> -e^ P ad u A pa , (43) 

the result for the axial current is more involved and includes the vector source 

_ 2 

1pl5l^ -> -e^padvBpv - — t^paSydpSv . (44) 

In our approach, this is a consequence of the anomalous behaviour of the 
fermionic measure under axial gauge transformations and in this way the 
bosonic form of the axial current correctly yields its anomalous divergence. 
We also mention the possibility of considering the particular case of a purely 
chiral external source [s^ = ±t^), and obtaining a bosonized version for 
this model. The Kalb-Ramond field then corresponds to a particular 'chiral' 
combination of A and B, namely Ap V ± B^ v . 



As stressed above, recipes (|43f) and (|4]) can be considered exact apart 
from the fact that if one is to work in the bosonic version one has to use an 
approximate expression for the bosonic action. The one we proposed is based 
in a quadratic approximation and leads to the bosonic generating functional 
presented in eqs. (^0D -(^2|). 



It should be noted that, besides playing an important role in the bosoniza- 
tion rule for the axial current, the chiral Jacobian also affects the actual form 
of the bosonized action, being the cause of the existence of non-quadratic 
terms in the currents, and of the coupling between the field strength for the 
Kalb-Ramond field B and the currents. This situation may be contrasted 
with the one of having just a vector current, where all the non-quadratic 
terms disappear if the fermionic determinant is evaluated up to second order 
in the fields, as we did. These non-quadratic terms are a signal of the anoma- 
lous Ward identity linking the divergence of the axial current and two vector 
currents. In spite of the fact that we haven't included the triangle diagram 
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in our approximate evaluation of the fermionic determinant, partial informa- 
tion from it has shown up from the Jacobian whose calculation implies the 
knowledge of the exact axial anomaly. It is interesting to note that, as it 
happens for the complete d = 2 bosonization recipe, the axial anomaly de- 
termines basic properties of d = 4 bosonization (For odd dimensional spaces 
it is the parity anomaly which seems to play a similar role 0-[0]). 

Finally, a comment about the choice of the actual value of the renor- 
malized 'mass term' for the axial source s M : If the source is not dynamical, 
there is no propagator associated to it and, needless to say, its natural renor- 
malized value is zero, since this value will not be modified by any other 
higher order correction (they would require diagrams with internal s M lines). 
When the source is dynamical, the actual value of the renormalized mass 
is arbitrary and becomes a new quantity to be 'measured'. A model with 
dynamics for s M is not necessarily anomalous, one may consider a system of 
many fermionic species, with their charges adjusted in order to cancel non- 
trivially the anomaly. It should be emphasized that if the vector source is 
not dynamical, the regularization can be chosen differently since the criterion 
of preserving the conservation of the vector current at the quantum level no 
longer applies. 
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